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Memory effects in transport require, for their incorporation into reaction diffusion investigations, a 
generalization of traditional equations. The well-known Fisher's equation, which combines diffusion 
with a logistic nonlinearity, is generalized to include memory effects and traveling wave solutions of 
the equation are found. Comparison is made with alternate generalization procedures. 
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I. INTRODUCTION 

Fisher's equation Q describes the dynamics of a field 
u(x,t) subject to diffusive transport and logistic growth: 



du ^d 2 u 



(1) 



This kind of reaction-diffusion equation is relevant in 
chemical kinetics as well as in ecological contexts where 
u diffuses via a diffusion constant D and grows with a 
linear growth rate k while the environment imposes a 
carrying capacity K . Equation (|lj) belongs to a family of 
single component models of broad applicability. Fisher 
proposed it as a deterministic model of the spread of a 
favored gene in a population . Other systems in which 
Eq. @) plays a significant role include flame propaga- 
tion, neutron flux in a nuclear reactor and the dynamics 
of defects in nematic liquid crystals |l], |[ |], ||, g] • Work 
has been done on the propagation of wave-like structures 
in unbounded systems jq, uj , as well as the formation of 
spatial structures in finite geometries ||. 

Diffusion is typically a limit of a more coherent mo- 
tion interrupted by scattering events which is valid when 
the scattering events are extremely frequent. A useful 
manner of describing situations in which the frequency 
of scattering events is not infinitely large is through the 
incorporation of memories with finite decay times H . An 
analytical investigation of reaction-diffusion phenomena 
based on the introduction of exponential memory in the 
context of the logistic nonlinearity of the Fisher's equa- 
tion, and a piecewise linearity approximation, was re- 
cently given by Manne, Hurd and Kenkre (MHK) ||o). 
A number of interesting results appeared in that analy- 
sis, including a generalization of the relation between the 
minimum speed of shock fronts and the system parame- 



ters, and the appearance of a speed limit beyond which 
the fronts would exhibit spatial oscillations. 

As stated in Ref. Efl|, the generalization of Fisher's 
equation is not unique, the one employed in Ref. fllOj hav- 
ing been chosen for analytic tractability. A more natural 
generalization results in a more complex structure of the 
equation, involves a nonlinearity in the damping, and is 
the subject of the present paper. Our investigations show 
that, in the light of the present generalization, some of 
the interesting predictions of Ref. |10| arc modified dras- 
tically, while others are only slightly changed. 

Reference pi first generalized the linear part of Eq. (fil) 
by incorporating an exponential memory, converted the 
resulting integro-differential equation via differentiation 
into a differential (specifically the telegrapher's) equa- 
tion, and finally added the nonlinear logistic term to ob- 
tain the starting point of the analysis. Here, however, we 
take as our point of departure, 



t-»jf*-'>5?*+*/M. 



(2) 



which reduces to Fisher's equation if the nonlinearity is 
logistic, and if the memory function <p(t) is a (S-function 
in time. Generally, the decay time of <fi(t) is a measure 
of the time between scattering events: we will take </>(£) 
to have a simple exponential form </>(£) — ae~ at , where 
1/a represents the scattering time. 

Transformation of Eq. (0) into a differential equation 
is trivially done by differentiating once with respect to 
time: 



d 2 u 
dt 2 



[a kf'(u)] ^ = v 2 ^ 



■akf(u), (3) 
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where, as in [ pD| , we put Da = v 2 , the physical mean- 
ing of v being the speed dictated by the medium in the 
absence of scattering. This is the speed at which the 
underlying quasiparticle, whose number density or prob- 
ability density is described by u(x, t), moves ballistically 
(coherently) in between scattering events. If the coherent 
motion is interrupted too often by scattering, the motion 
looks diffusive and one returns to the Fisher limit. 



Equation (||) becomes, in the absence of the nonlin- 
earity, k = 0, the well-known telegrapher's equation |ll]] 
suggested by Lord Kelvin for the description of propa- 
gation of transatlantic telegraph signals. Equation (J3|) 
differs from the starting point of the MHK analysis, viz., 



8 2 u 
~dH 2 



du 



d 2 u 
dx 2 



akf(u), 



(4) 



in a significant way. The former shows that the damping 
coefficient multiplying du/dt is not constant, and can 
even be negative. 

Whereas Eq. (|f|) is obviously simpler in form than 
Eq. (0), the latter could be argued to stem from a more 
natural generalization of the Fisher equation, involving a 
simple replacement of a S- function memory by a finite de- 
cay time memory in the transport term Dd 2 u/dx 2 in (Q). 
We are interested in comparing the physical consequences 
of the two generalizations in the context of traveling wave 
solutions. We do this in two steps. In Section O, we give 
a general analysis of traveling wave solutions without ap- 
proximating the equation in any way. We investigate the 
minimum speed of wavefronts, the existence of spatially 
oscillatory fronts given in MHK, and give a qualitative 
description of all the kinds of front shapes that the sys- 
tem can support as a function of the parameters. In a 
certain regime in parameter space, a dynamical stabi- 
lization of the unstable state of the logistic equation is 
observed, with fronts of state u = invading the state 
u = K. In Section III, we replace the nonlinearity by 



a piecewise linear approximate form and obtain analytic 
solutions following the method of MHK. We end with 
concluding remarks in Section IV. 



write it as a first order system: 

U' =V =f(U,V) 

V =mk-a-2kU)cV-akU(l-U)} = g(U,V). 

(7) 
The system (0) has two equilibria: (U*,V*) = (0, 0) and 
([/*, V*) = (1, 0). We can analyze the character of these 
by looking at the linear behavior in a neighborhood of 
the equilibria. 

At the equilibrium (0,0) we have the eigenvalues: 



A, 



c(k - a) =F \/(k - a) 2 c 2 - Akc 
2m 



(8) 



If we assume that U is a density or a concentration, so- 
lutions where U oscillates below U — are not allowed. 
This imposes a condition on the eigenvalues Azp to be 
real, from which we obtain: 



C -~* Cmin. ^ 



f 



l + My-1/yf 



(9) 



where y = ya/k. This relation states that there is a 
minimum value of the speed of the nonlinear waves. It 
is to be compared to the known result in the context of 
Fisher's equation: c m in = 2\fkD (see for example M). If 



we make the formal identification Da — v , as in (H) , we 
have: 



c > c„ 



(10) 



Equation (0) is also to be compared to the previously 
found result of Rcf. [Hoi: 



II. NONLINEAR ANALYSIS 

We look for waves moving in the direction of increasing 
x: u(x, t) — K U{x — ct) = K U{z), where c is the speed 
of the nonlinear wave, generally different from the speed 
of linear waves v, dictated by the medium. We are also 
renormalizing with respect to the carrying capacity K 
to simplify the reaction term in the logistic case. With 
this ansatz, we obtain from (H) the following ordinary 
differential equation: 

(v 2 - c 2 ) U" +c[a- kf(U)} U' + akf(U) = 0. (5) 

For a logistic reaction term: /(£/) = U(\ — U) and 
f'(U) = 1-2U, so that Eq. (|) becomes: 

mil" + c(a-k + 2kU) U' + ak 17(1 -U) = 0, (6) 

where we have followed the notation of MHK, 771 = 1; — 
c 2 , to emphasize the formal similarity between (|6|) and 
the equation of motion of a damped oscillator of mass m, 
subject to the nonlinear force — akU(l — U). 

A convenient way to analyze the solutions of (pi) is to 



> MHK _ 



1 



1 



\v 2 



(11) 



Figure |l| displays a comparison of the minimum speeds 
Cmin as a function of the system parameter a/k in the 
three cases: purely diffusive (Fisher limit), the MHK gen- 
eralization, and the present generalization. It can be seen 
that the two generalizations approach asymptotically the 
behavior of the purely diffusive situation for large values 
of a/k. This is to be expected since a — > 00, v — > 00, 
v 2 /a = D is the diffusive limit of Eq. (0)]. At low values 
of the ratio a/k, however, there is a notable difference. 
The finite correlation time generalizations allow waves 
with lower speed than those present in Fisher's equation. 
Moreover, the present generalization predicts a sharply 
different behavior in the region a < k, due to the anti- 
damping present in Eq. (|6|). This growing branch of the 
c-min function turns out to have drastic consequences in 
the nature of the traveling fronts when the system is far 
from the diffusive regime. On the one hand, it can be 
seen that waves with speed c = v are obtained at a finite 
value of the memory constant a, namely a = k. When 
approaching the ballistic limit beyond this point, trav- 
eling waves are allowed with speeds lower than v. The 



1.0 



with memory (this generalization)" 
with memory (MHK) 
purely diffusive 




FIG. 1: The ratio of the minimum speed of nonlinear waves 
Cmin to the medium-dictated speed v, as a function of the 
system parameter a/k, which is the ratio of the scattering 
rate at which motion coherence is interrupted to the growth 
rate. Lower values of a/k represent more coherent motion. 
Three curves are shown: the prese nt g eneralization (full line), 
the MHK generalization of Ref. [ncl (dashed line), and the 
purely diffusive Fisher limit (dotted line) with D = v 2 /a. 



nature of these will be analyzed in the following discus- 
sion. 

Around the equilibrium (1,0) the eigenvalues are: 



Mt 



-c(k + a) =F y(k + a) 2 c 2 + Akam 
2^ ' 



(12) 



Oscillations around this equilibrium are in principle pos- 
sible, as observed in 1 10 , since there is no impediment to 
the solution growing beyond U = 1. The condition for 
oscillations is found from the radicand of the eigenvalues 
/i T : 



(k + a) 2 c 2 + 4ka{v 2 - c 2 ) < 0, 
which, after some manipulation, becomes 
(a - k) 2 c 2 + Akav 2 < 0. 



(13) 



(14) 



Both terms in this relation are positive, implying that the 
relation cannot be satisfied by the wave speed for any set 
of parameters. 

A similar analysis can be carried out with the equa- 
tion ([|), corresponding to MHK model. The character 
of the equilibria for both generalizations is resumed in 
Tables | and O. Note that there is a c osc in Table || Case 
(7) , not present in the results resumed in Table B because 
Eq. ( |l4| ) cannot be satisfied. Different combinations of 
parameters allow a variety of traveling fronts connecting 
one equilibrium to the other. Cases (1) to (7) in Table || 
are the fully nonlinear equivalents of the piecewise linear 
situations analyzed in [fill. The situations resumed in 
Table R are discussed below. 



Case No. 


a, k 


c 


(0,0) 


(1,0) 


(1) 


a > k 


C *C Cmin 


stable spiral 


saddle 


(2) 


Cmin \ C <^ V 


stable node 


saddle 


(3) 


V < C 


saddle 


unstable node 


(4) 


a < k 


C *C Cmin 


unstable spiral 


saddle 


(5) 


Cmin \ C <i, V 


unstable node 


saddle 


(6) 


V < c 


saddle 


unstable node 



TABLE I: Character of the equilibria for the possible com- 
binations of parameters and wave velocity, according to the 
generalization studied in this paper. The "case number" refers 
to the discussion in the text. 



Case No. 


a, k 


c 


(0,0) 


(1,0) 


(1) 


a > k 


C \ Cmin 


stable spiral 


saddle 


(2) 


Cmin *C C <C V 


stable node 


saddle 


(3) 


V < C 


saddle 


unstable node 


(4) 


a < k 


C \ Cmin 


stable spiral 


saddle 


(5) 


Cmin \ C \ 1 


stable node 


saddle 


(6) 


V < C < C 03C 


saddle 


unstable node 


(7) 


Cosc V, C 


saddle 


unstable spiral 



TABLE II: Character of the equilibria for the possible com- 
binations of parameters and wave velocity, according to the 
generalization studied in jlOJ. 



Equation. (|6|) can be exploited directly to understand 
the nature of the solutions of the traveling wave prob- 
lem, without obtaining analytical solutions via simpli- 
fications. We will base the following discussion on the 
mechanical interpretation of Eq. ml), which can be rein- 
terpreted, following MHK, as describing the motion of a 
particle of mass to in a nonlinear potential, subject to a 
state-dependent damping that can be positive or nega- 
tive. Since the mass of the particle can also be positive 
or negative, depending on whether the velocity c is lower 
than v or not, we analyze the two cases separately. 



A. c < v 

In this case a particle of mass m > is moving in 
a potential <f>(U) = ak(U 2 /2 - U 3 /3) [see Fig. § (a)]. 
This potential has a minimum at U = 0, but we have 
ruled out solutions that oscillate around it based on the 
positivity of U [this are cases (1) and (4) in Table pi. 
Another possible solution is an overdamped trajectory 
that connects the equilibrium at U — 1 with that at 
U = 0. This corresponds to a traveling front of the state 
1 invading the system at state 0. The damping coefficient 
is 7 = c(a — k + 2kU). We can see that if a > k, then 
7 > for all values of U along this trajectory. This is 
a solution that connects the saddle with the stable node 
shown in the case (2) of Table |[ 

Now, if a < k, the motion is damped on part of the 



trajectory [when U > (1 — a/k)/2, see Fig. g (c)] and 
anti-damped on the rest, specifically near [7 = 0. This 
means that the equilibrium at U = is unstable, and 
there is a solution connecting it to the equilibrium at 
U = 1. It is a trajectory connecting the unstable node 
at (0, 0) to the stable manifold of the saddle at (1, 0) [see 
Table B, case (5)]. This trajectory corresponds to a front 
of state invading the state 1. 



B. 



v < c 



In this situation m is negative, so we multiply Eq. 
by — 1, to obtain: 

\m\U" -c{a-k + 2kU)U' = akU(l - U). 



(15) 



Now the particle of mass \m\ moves in the potential 
</)(U) = -ak(U 2 /2 - U 3 /3) [refer to Fig. § (b)], where 
there is a stable equilibrium at U = 1. The sign of the 
damping coefficient is also reversed, and the motion is 
always anti-damped when a > k. When a < k, it is 
always anti-damped in the vicinity of U = 1. It can be 
seen in Table | [cases (3) and (6)] that the equilibrium 
at (0, 0) is always a saddle, and the equilibrium at (1, 0) 
is always an unstable node. There is a trajectory that 
connects the unstable node to the stable manifold of the 
saddle, corresponding to a front of state 1 invading the 
state 0. 

We note that, contrary to the generalized telegrapher's 
equation studied in pffl, we do not have here oscillat- 
ing wave shapes (corresponding to the unstable spiral at 
(1,0) shown in Table |J). Also to be remarked are the 
solutions where the state U = — which is an unstable 
state of the logistic equation — invades the state U = 1. 
This stabilization of the null state has been made pos- 
sible by interplay of memory and reaction, through the 
state-dependent damping coefficient of Eq. (0) . 

A phase diagram summarizing these results is shown 
in Fig. |3| A qualitative indication of the nature of the 
front wave is given, together with a reference to the cases 
enumerated in Table \ 




FIG. 2: The potentials and the damping coefficients of the 
mechanical analog of Eq. (h). Cases (a) and (c) correspond 
to m > 0, while (b) and (d) correspond to m < 0. 



responds to a = b), the oscillator equations in the trav- 
eling waves ansatz (|5j) become: 

mU" + 2 lx U' + k\U = 0, U<a, (17) 
to U" + 2 72 U' + k\ {b - U) = 0, U>a. (18) 



where 



c / k\ c 

7i = 77 " , 72 = ^1" 



b — a I 



k\ 



ak 
a 



ak 



(19) 



(20) 



Suppose that we are looking for a front solution that 
interpolates from U = b as z — ► — oo to U — as z — > oo. 
Without loss of generality, we take U — a at z = 0. 
We look for solutions of two oscillators: one to the right 
of z = 0, and one to the left of z = 0. We define the 
following variables: 



III. PIECEWISE LINEARIZATION 

Explicit solutions corresponding to the front waves de- 
scribed above cannot be found in the fully nonlinear sit- 
uation. However, a piecewise linearization of Eq. (||) can 
be made, just as in 1 10 , taking the following as a reaction 
function: 



f(U) 



U/a, U <a, 

(b-U)/(b-a), U>a, 



(16) 



where a < b. With this reaction term, which inciden- 
tally, generalizes somewhat the logistic form (which cor- 



U R (z) = U(z), z>0, 


(21) 


U L (-z) = b-U(z), z>0. 


v" J v 


7|) and (|18|) become: 




to U& + 2 7l U' R + k\ U R = 0, 


(22) 


-mU'l + 2 72 U' L + kfU L = 0. 


(23) 



Equations ([1 



The wave front shape, which is a solution to Eqs. (O) 
and (|l8|), can be easily found in terms of exponentials, 
by matching solutions of Eqs. (22) and (23) and their 
derivatives on either side of z = 0. In some cases, one 
of the solutions of Eqs. p^) or ( p3[ ) contains a grow- 
ing exponential, and the corresponding factor has to be 
set to zero to avoid unbounded growth. Examples of all 



IV. CONCLUSIONS 




FIG. 3: The regions in parameter space where different trav- 
eling fronts are found. The front shapes are schematic. The 
numbers refer to the cases shown in Table |. 



U 0.4 




FIG. 4: Front waves in the piecewise linear version of the 
model, a < 2k and a > 2k correspond to the two distinct 
regimes of damping or anti-damping, in this linearized model. 
We have taken a = 1/2 and 6=1. 



the cases described in the previous section are shown in 
Fig. El within this piecewise linear scheme. The solutions 
have been arbitrarily shifted in the variable z — x — ct 
to avoid superposition of the curves. In the model lin- 
earized according to (|16| ), the regimes of damping and 
anti-damping of Eq. rtq) are separated by the condition 
a = 2k, which is used in the figure caption to classify 
the curves. (The condition separating the damping from 
the anti-damping regimes is a = k/a in the piecewise 
linear model, instead of the a — k of the fully nonlinear 
model.) The case numbers corresponding to Table Q are 
also shown in the caption. 



Reaction-diffusion systems in which the transport pro- 
cess is wave-like at short times and diffusive at long times 
are the focus of the present investigation. This passage 
of the character of the motion from coherent to incoher- 
ent is a general feature of all physical systems and may 
be represented by a memory function whose decay time 
(or correlation time) represents the demarcation. Wc 
have analyzed here a generalization of Fisher's reaction- 
diffusion equation which has a logistic nonlinearity de- 
scribing the reaction process. Our study has centered on 
traveling wave solutions. From arguments without ap- 
proximation we have found a generalization of the known 
Fisher's equation result regarding the minimum speed of 
the traveling waves. While the generalization converges 
to the Fisher's equation result in the limit of diffusive 
transport (memory that decays infinitely fast), sharp dif- 
ferences occur in the wave- limit (see Fig. [fl): new kinds 
of solutions are found to be possible involving "inverse" 
fronts, in which the state [7 = invades the state U = 1 
(see Fig. [|). We have also found explicit analytic solu- 
tions (typical cases plotted in Fig. [|) via the piecewise 
linearization introduced for this problem elsewhere |1C[| 
and pointed out a number of differences in the predictions 
arising in our present analysis. Some of these differences 
are expected but at least one of them is surprising: so- 
lutions which oscillate spatially when they have a speed 
above a certain value (called c osc in Ref . pi ; see also Ta- 
ble |n|) predicted in Ref. M] are found to disappear in the 
present analysis which can be argued to be a more nat- 
ural generalization of Fisher's equation to include wave- 
like transport. There are a number of issues, such as 
the stability of the solutions against perturbations intro- 
duced into the system, which will be studied in future 
investigations. 

We mention in passing the problem of what kind of 
initial conditions eventually evolve into a traveling front. 
The general problem can be very difficult, but a simpli- 
fied treatment is possible as follows. Consider that the 
initial state has the behavior u(x, 0) ~ Ae~ ax for x — > oo. 
Correspondingly, we suppose that the leading edge of the 
traveling wave has the form: 



i{x,t) = Ae- a{x ~ ct \ 



(24) 



with a and A arbitrary. 

We substitute this into the differential equation, and 
suppose that, at the leading edge, u s=s 0: 

(ac) u + aac u — kac u + o(u ) = v a u + ak u + o(u ) . 

(25) 
Disregarding the terms in u 2 we arrive at a dispersion 
relation between c and a: 



2a 



n j(a + k) 2 



Aa 2 



(26) 



The negative sign before the square root in this expres- 
sion leads to negative c, and can be disregarded. The 



expression with the positive sign gives the speed of the 
front wave that will eventually develop from an initial 
condition that has an exponential decay of coefficient a. 
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